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ABSTRACT 

In this paper, we have applied a new Chebyshev wavelet based coupled method to film-pore diffusion model. 
We have used the film-pore diffusion model, satisfactorily describe kinetics of methylene blue adsorption onto the three 
low-cost adsorbents, Gauva, teak and gulmohar plant leaf powders. We first discretize the nonlinear PDE and then 
utilized the Chebyshev wavelets on the obtained discretize nonlinear differential, equation to get the solution at each 
approximation of Picard method. The proposed wavelet based results are compared with the results, obtained by method 
of lines (MOL) and experimental results. Satisfactory agreement with the experimental and MOL results is observed. 
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1. INTRODUCTION 

Large amounts of dyes are utilized by textile industry and a significant portion of these dyes are not 
consumed in the process and therefore let out with the effluent. In recent years, the cost of commercial adsorbents 
is too high; interest for using low-cost adsorbents for removal of dyes from textile effluents is continuously 
growing. A recent survey indicates that, in India, on an average fresh water consumed and effluent generated per 
kg of finished textile are 175 L and 125 L, respectively [1] [Ponnusami et al 2010]. The presence of dyes in 
aqueous effluents is highly objectionable, as this affects the photosynthetic activity in receiving water body by 
reducing/preventing light penetration. As the dyes are recalcitrant in nature it is difficult to treat them, in 
conventional biological treatment plant [2, 3] [Ponnusami et al. 2007, Ponnusami et al. 2009]. Various researchers 
have worked on biological degradation of dyes. But, very often, the metabolic intermediates are found to be more 
toxic than the original compound [4]. Therefore, identification of low-cost adsorbents is given more attention by 
the researchers recently, as commercial adsorbents like activated carbon are too costly. Few recent studies 
investigating application of low cost adsorbents are: jackfruit peel [5], pine apple stem [6], phoenix tree leaves [7], 
pomelo peel [8], shells of bittim [9], orange peel [10], broad been peels [11] etc. 

In our previous research works, we investigated the feasibility and adsorption of MB onto three plant leaf 
powders namely guava leaf powder (GLP), teak leaf powder (TLP) and gulmohar leaf powder (GUL) [1]. 
Hariharan et al. [46], used the film-pore diffusion model (FPDM) to describe the kinetics of methylene blue 
adsorption onto GLP, TLP and GUL. Diffusion based kinetic models are too complex and they require rigorous 
solution methods. For many of the diffusion and kinetics models do not have analytical solutions. In our previous 
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work we had used method of lines to solve film-pore diffusion model and had shown that Film-pore model could describe 
the kinetics of adsorption of MB onto GLP, TLP and GUL [1], Hariharan et al. [46], had applied the Haar wavelet method 
for the solution of the above model. In this work, we have proposed a wavelet based coupled method to film-pore diffusion 
model. 

Recently, wavelet based methods become important tools, because of the properties of localization. Chebyshev 
wavelet method (CWM) is one of the earliest and simplest examples of what is known now as a compact, dyadic, 
orthonormal wavelet transform [12-24,26]. Recently, Hariharan and his co-workers [19-32] have applied the wavelet 
transform methods for solving differential equations arising in Engineering. Wavelet-Galerkin method (WGM) was 
successfully applied for solving the nonlinear differential equations [33, 34]. 

The basic idea of Chebyshev wavelet method is to convert the PDEs to a system of algebraic equations by the 
operational matrices of integral or derivative [35-45]. The main goal is to show how wavelets and multi-resolution analysis 
can be applied for improving the method in terms of easy implementability and achieving the rapidity of its convergence. 
In this work, we have applied a wavelet based coupled method, which combines the Picard’s method and the Chebyshev 
wavelet method for the numerical solution of the film pore diffusion mode. To the best of our knowledge until now there is 
no rigorous Chebyshev-Picard wavelet solution has been reported for the above model. 

2. MATERIALS AND METHODS 

Film-pore diffusion model (FPDM) is successfully described earlier by McKay and co-workers [47, 48]. 
Numerical results of FPDM by method of lines are described in our previous paper [1]. Haar wavelet based solution of 
FPDM is described in our previous paper [46]. In the present paper, we have developed the Chebyshev-Picard wavelet 
method is described in detailed and the results are compared with our previous solution. 

2.1. Wavelet Preliminaries 

Wavelets are the family of functions which are derived from the family of scaling function { 0j,k-. k Z J 

where: 

0(:r) = a k 0( 2x — k} (2.1) 

For the continuous wavelets, the following equation can be represented: 

(x) = (~) £ i?, a =£ On (2.2) 

Where a and b are dilation and translation parameters, respectively, such that ¥ f:r) is a single wavelet function. 
The discrete values are put for a and b in the initial form of the continuous wavelets, i.e.: 
a = a 0 "- i , a 0 > 1, b 0 > 1, (2.3) 

b=kb Q a Q * r j,k GZ. (2.4) 

Then, a family of discrete wavelets can be constructed as follows: 
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Yjik = x ~ *)' ( 2 . 5 ) 

So, v constitutes an orthonormal basis in L 2 (R), where (P (jf) is a single function. 

2.2. Function Approximations and Chebyshev Wavelets 

We can express any function f (A) £ Ls [(), 1) into truncated second kind Chebyshev wavelets series as 
indicated below. 


2 l ~' M -1 

A function / (*) = £ j] C m \f/ nm (x) = C Tx ±> (x), 

n =1 m—0 


( 2 . 6 ) 


where C and v P(jc)are 2* *Mxl matrices, given by the following column matrices. 


c 




(P 10 [x),y/ n (x),...,^6 M _, (x),^ 20 (x),...,^ 2jW _ 1 (x),...,y/ 2k o (f),...,^ 2 i M1 (a-) 


Using the collocation points as following: 




2 fc M 


M. 


We can define the Chebyshev wavelet matrix V P H 


as: 


*P □ 

fflXm 


*P 


YM 


, 'P 


YM 


.... *P 


^ 2*M-1 ^ 

v YM j 


For example, when M = 3 and k = 2 the Chebyshev wavelet is expressed as 


1.2732 

1.2732 

1.2732 

0 

0 

0 

0 

0 

0 

1.2732 

1.2732 

1.2732 

-1.6977 

0 

1.6977 

0 

0 

0 

0 

0 

0 

-1.6977 

0 

1.6977 

0.9903 

-1.2732 

0.9903 

0 

0 

0 

0 

0 

0 

0.9903 

-1.2732 

0.9903 


2.3. Chebyshev Wavelets Operational Matrix of Integration 

The integration of the vector *P (f) defined in Eq. (8) can be obtained as 
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J l F(s)dsD PW(t), 

0 

where P is the ( 2 k M \ X ( 2 k M ) operational matrix for integration and is given 


(2.7) 


as 


f c s s 
o c s 
o o c 


s^ 

s 

s 

s 

c 


0 0 0 

V / 

where S and C are M/M matrices given by: 


5 = 


and 


f 

1 

0 
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o' 


0 

0 
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1 

0 
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~3 






0 

0 
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1 

0 
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1 

0 
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^ 2M(M-2) 






] 
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2-Jl 

0 
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0 

1 
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8n/2 
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1 

8 

1 

C= i 

i 

&/2 


1 

~4 

0 

1 

12 

2 







2j2(M-l)(M-3) 

1 


4(M-3) 

0 


4(M-2) 


4(M-1) 

0 


The integration of the product of two Chebyshev wavelet function vectors is obtained as 


= j v F(.T) v i /r (x)J.x 


( 2 . 8 ) 


where I is an identity matrix. 

Theorem 2.1 (Accuracy estimation): Let /(f) be a continuous function defined on[0, l), with bounded 
second derivative |_/"^(T )| bounded by B , and then we have the following accuracy estimation: 


^k,M < ' 


4kb 


y -Ly_ i — 

/ j 5 / 4 / \ 4 

«=/+ i n m~M ym 1 j 


(2.9) 
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where 


<7 k.M 


p M -1 


J / (') “ Z Z C n m ¥ ( n,l {<) W n (t)clt 


n =1 m =0 




2.4. Convergence Analysis of the Proposed Method 

In this section, we can derive an error bound of the Chebyshev wavelet Picard technique to an arbitrary unknown 

function. 


Theorem 4.1: Let S, k and M tends to °°, then the series solution Equation.(9) converges to f (A') . 
Proof: Let L~ [0,l) be the Hilbert space and y/ n m forms a basis of L~ ! [o,i) . Consider 


2 W M- 1 

/ I+ i(- i ) = XZcv,,„(4 

n =1 ni—0 

where 

l 

Cnt = {fs + 1 M . Wn.m (*)) = / f (*) W n , n (*) dx, 

0 

Let R s k + h be a sequence of partial sums of C s ^\j/ n m (t) , we prove that Rj is a Cauchy sequence in Hilbert 
space L" [0, l) and then we show that Rf. + h converges to f s+1 (%) , when k. M -A 

We show that is a Cauchy sequence. Consider R' k + h be arbitrary sums of with 

k > k, M > M. 


1 _ r>v+l 


2 M-l 


2 M-l 


n=l m=0 
2 


Z Z (*) - Z Z (*) 

1 m—0 n= 

2 t-1 M-l 

Z 

:2 t_1 +1 m-M 
2 k ~ l M-l 

= ( z Z C ™V„,„(^)> z Z c rVu( x ) 

\«=2 W +1 m=M i=2 k ~ l +1 j=M 

2 k ~ l M-l 2 l_1 M-l 

= z z z Z c » c r(^w>^w 


M-l 


n=2 k ~' + \ m-M i=2 k ~'+l j=M 
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M-l 


E El 
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Using the Bessel’s inequality, we have 


ZZb 

n—l m =0 


2 


is convergent and hence 


nJ+l _ pi+1 


■0 as k,M, k,M —> °° . This implies that R, y is a Cauchy sequence and it converges to 


T i+1 (.x:)e L 2 [0,l). It is enough to show that y j+1 (.*) = f s+l (.*) . 


5+1 5+1 

= r — r 

nm nm ’ 

=0 


2 1-1 M-l 

Hence ZZ c «mV«.m( x ) converges to / s+1 (.*) as 

n—l m—0 


Hence the theorem 

3. METHOD OF SOLUTION BY THE CHEBYSHEY WAVELETS- PICARD METHOD (CWPM) 

The film pore diffusion model is governed by a Reaction-diffusion-convection equation [1,46] 

dc A <^\d 2 C ldc 

— = A(C)—+ -— (3.1) 

dr dz~ z oz 

Where 


A(C) = 


£ P + 


(QhP P 

f . , \ 

1 + bc 0 

V C 0 y 

v (l + bc 0 Cf p 


with initial conditions 

C(z,0) = , C(z,l) = e~ z - 0 09 , C(0, t) = e~ 0X>9r , c( 1, t) = e" 1 " 0 097 (3.2) 


C -Concentration of the dye (methyline blue), T -time variable, z-spatial variable, £ p -particle porosity, p - 
particle density C 0 — initial bulk concentration 

We apply the Picard technique to equation Equation.(3.1) 


dC 


5+1 


-A(C) 3=C "' 


iSC, 


Sr Sz 2 z Sz 

Where SG N ,C (z,T) is known and this term can be used for obtaining C +1 (z,T) 


(3.3) 


The initial and boundary conditions are 
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„- 1 — 0.097 


C J+1 (z,0) = e-\C s+1 {z,l) = e- z -° 09 ,C s+1 {0,T) = e^ m \C s+l {\,r) = 

We can establish the Chebyshev wavelets method (CWM) to equation Equation.(3.3) 

m m 

= S S CljY l (z)yAz)= y T (z)F s+ 'y/{r) 


d 2 C s+x 


dz 

Apply the integral operator In 7 Equation.(3.4) 

C J+1 {z, t) = In:y/ T (z)F i+ V(r) + /,. (t)z + 1 2 (r) 

Using the boundary conditions, 

1 2 (t)= 0 

b 

Equation.(3.5) becomes 

c s+1 {z,r) = ln:y/ T {z)F s+l ys{r)-x[lni h ys T {z)]f S+I ys{r) 

Equation.(3.3) can be written as 


dC. 


= A(C) ( - 


d 2 C s+l 1 dC s 


+ 


dr dz 2 z dz 

With the help of Equation.(3.4), Equation.(3.7) becomes 


dc 


5+1 


dr 


■ A(c)i/s T {z)F s+x y/{r)-y/ T (z)iV' 5+ y(r) 


where AT 1 = (z)) ‘ 


(3.4) 


(3.5) 


(3.6) 


(3.7) 


(3.8) 


(3.9) 


Using the integral operator In L to Equation.(3.9) and the initial conditions to obtain the following formula 

c s+l {z, t) = A(c)y/ T {z)f s+x (ln' L y/(j))-y/ T ( z)N’ +1 ^ + e“ z “ 009 

From Equation.(3.6) and Equation.(3.10), one can get 

G[ln 2 z¥ T {z)-z[lnl b y T {z)]F s+1 - F s+l \l l L i/s (t)(i/s(t )) =1 + GFys{r)]' J= 0 
Where, G = (a(c)i/s t (z))“‘ and F = [- (- y T (z)(V' +I (l l L i/T{r))+e- z + e“ z " 009 )] 


(3.10) 


(3.11) 


Equation.(3.11) is the Sylvester equation. We solve Eq.(3.11) for F s+i and applying F'' +l in equations (3.6) 
(3.10) to obtain the solution C{z, 7") at (.S’ +1) approximation of Picard iterative approach, at the collocation points. 
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Table.1: Comparison of CWPM, MOL and HWM by Obtaining the Mass Transfer Coefficients 
Using Film-Pore Diffusion Model Adsorption of MB onto GLP and k=l and 
M=3. [ A(Cj) are Constants (linear) and e = 0.5, p = 500] 


Temp 

c ° 

( mgdnt ) 

External- film transfer coefficient 

A ) 

Internal effective diffushitv 

Error 

K 


MOL 

HWM 

CWPM 

MOL 

HWM 

CWPM 

Em 

E, 

Ec 

303 

50 

100 

150 

200 

1.00x10* 

3.00x10* 

4.23x10* 

1.74x10"“ 

2.54x10"“ 

1.24x10“ 

1.197 

0.140 

0.935 

1.61 

0.938 

0.129 

0.824 

1.102 

0.638 

0.029 

0.483 

0.873 

313 

50 

100 

150 

200 

1.71x10* 

3.45x10* 

6.35x10* 

6.46x10"“ 

9.26x10"“ 

7.32x10“ 

1.462 

1.120 

1.267 

7.570 

1.221 

0.927 

1.016 

5.112 

1.112 

0.653 

0.820 

3.234 

323 

50 

100 

150 

200 

4.27x10* 

4,75x10* 

4.01x10* 

3.11x10"“ 

5.3Ixl0 13 

4.54x10“ 

0.856 

0.160 

2.168 

3.164 

0.28 

0.000 

1.208 

1.016 

0.192 

0.000 

0.727 

0.523 


4. CONCLUSIONS 

In the present paper, we have successfully applied the CWPM for the film-pore diffusion model. It was 
investigated that the model could predict the concentration decay curve for all adsorption of methylene blue onto TLP, 
GUL and GLP excellently with a small deviation, during initial period. It gives better results when compared with the 
results by method of lines and Haar wavelet method. It is worth mentioning that, CWPM provides excellent results, even 
for small values of m . For larger values of ill , we can obtain the results closer to the real values. The method with far less 
degrees of freedom and with smaller CPU time, provides better solutions than classical ones. The work confirmed the 
power of the CWPM, in handling nonlinear equations. This method can be easily extended to find the solution of all other 
non-linear diffusion equations too. 
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